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Letter to the Editor 

Evidence for an inner molecular disk around massive Young 

Stellar Objects * 
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Abstract. We present observations of CO overtone bandhead emission toward four massive Young Stellar Objects 
(spectral type 06-B5). The high signal-to- noise ratio if-band spectra were obtained with VLT-ISAAC at a 
resolution Av — 30 km s _1 , sufficient to resolve the bandheads, but not the individual J- lines. We are able to 
explain the shape of the lines by assuming a simple isothermal keplerian disk model seen at different inclinations. 
The gas temperature ranges from 1500 to 4500 K and the CO column density is between 0.1 and 4 x 10 21 
cm -2 . The emission probably arises within the first few astronomical units of the disk, consistent with the high 
gas temperature. Our results indicate that molecules can survive close to a hot star and suggest that dense 
(n H > 10 10 cm" 3 ) inner disks may be relatively common at an advanced stage of high-mass star formation. 
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1. Introduction 

Disks are found around many young low-mass stars as 
a natural byproduct of angular momentum conservation 
during the star formation process. Observational evidence 
for disks around their more massive counterparts remains 
sparse, although they may be essential d uring the forma- 
tion of high-mass stars l|Churchwell2002|) . Extended disks 
have been found around massive young stellar objects 
(YSOs) in the Hot Core phase (e.g. IShepherd et alJl200li 
Beltran et all 120041: iGarav fc Lizanol Il999t IChini et all 
20041) . albeit those millimeter observations solely probe 



gas beyond 1000 AU. In the later Herbig Be stage, dis ks 
have only been recently discovered l|Fuente et al.l l2003h 

Evidence for the presence of a gaseous disk around 
the early B star NGC2024-IRS2 is presented by 
iLenorzer et alJ ((2004) . They argue that the infrared ex- 
cess observed toward NGC2024-IRS2 is best explained by 
the emission from a dense gaseous disk as close as 0.3 AU 
to the star. 

At the high densities and temperatures (1000-2000 K) 
characteristic of disks at 0.1-5 AU around young stars, 
molecules are expected to be sufficiently excited to pro- 
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* Based on observations collected at the European Southern 
Observatory at La Silla and Paranal, Chile (ESO Programmes 
68.C-0652 and 69.C-0448) 



duce a rich ro-vibrational spectrum in the near- (over- 
tone transitions Av — ±2) and mid-infrared (fundamen- 
tal transitions Aw = ±1). CO first overtone emission is 
most likely emitted by a disk an d /or a wind and ha s 
been detected and a nalyze d by, e. g . JScoville et alJ lTl979lh 
iGeballe & Perssonl il98<t : ICarrl lll989lk IChandler et alJ 
lll993lll995h" toward a number of other high- and low-mass 
stars. Their studies were, however, hampered by the small 
number of observed J-lines or a low resolving power and 
mo derate signal-to- noise . 

iBik et ah! l)2004[) examine -RT-band spectra of a sample 
of massive YSOs associated with ultra-compact Hn re- 
gions (UCHII) . Here, we present a simple disk model that 
describes the CO bandhead emission detected in four of 
these objects. The observations are presented in Sect. [21 
The disk model and the results are shown in Sect. [21 In 
Sect. 0] we discuss the location and origin of the CO. The 
implications of our findings on the process of high-mass 
star formation are summarized in Sect. [31 

2. Observations 

High signal-to-noise ratio long-slit i^-band spectra of four 
massive YSOs have been obtained in 2002 with the ISAAC 
spectrometer mounted on the VLT. A narrow slit (0.3") 
is used, resulting in a resolving power of R = 10 000 (Av 
= 30 km s _1 ). The spectra were reduced using standard 
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Fig. 1. If -band observed (upper spectra) and modeled (lower spectra) spectra of the CO 2-0 bandhead. The name of 
the object and the derived inclination is displayed. By definition 0° inclination corresponds to a disk seen pole-on. 



procedures. The wavelength calibration is based on the 
telluric OH emission line spectra. Telluric standard stars 
of spectral type AOV, observed under identical sky condi- 
tions as the science targets, were used to correct for the 
telluric absorption lines. A detailed desc ription and analy- 
sis of the if -band spectra can be found in lBik et~al . (2004). 
The properties of the objects are summarized in Tabled 

3. Disk model 

Continuum subtracted and normalized first-overtone 
bandhead spectra are shown in Fig. ^ (upper spectra). 
Four overtones are detected but only the first one is dis- 
played. The shape of the bandhead varies from source 
to source. The rise of the first bandhead is sharp in the 
spectrum of 08576nr292 and 11097nrl218 whereas in that 
of 16164nr3636 and 18006nr766 a blue wing is clearly 
present. 

The expected profile from a spherical stellar wind has a 
flat-topped shape and cannot account for the wing seen in 
two of the objects. The favored model to explain the blue 
wing consists of a keplerian disk and/or a disk- wind (e.g. 



IChandler et"ai]ll99R|) . although the disk model provided 
the best match. In this paper, we explore the possibility 
that the CO bandhead emission arises solely from a small 
disk in keplerian rotation around the massive star. We 
generate synthetic CO bandhead spectra using a standard 
parametric disk model for a range of gas temperatures, 
column densities, turbulen t velocities and disk viewing an- 
gles (cf. iKraus et alJl2000l for a detailed description). In 
contrast to standard disk models, the disk is assumed to 
be entirely isothermal. The molecular data (line frequen- 
ci es and Einstein c oeffici ents) are taken from the database 
of IChandra et al.1 l|l996l) . At temperatures above 1500 K, 
needed to populate the higher vibrational levels, no dust 
grains can survive and, therefore, no dust opacity is in- 
cluded. The population of the CO rotational levels within 
each vibrational level is assumed to be in local thermo- 
dynamic equilibrium (LTE), at the same temperature T ex 
as the local vibrational temperature. Each bandhead is 
characterized by a specific vibrational temperature. The 
downhill simplex method is used to determine the model 
parameters: the CO total column density A^(CO), the ex- 
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Fig. 2. Fit of the complete observed CO bandhead spectrum of 11097nrl218. The upper curve is the observed spectrum 
whereas the lower one is the model spectrum. 



IRAS 


Object 


Mass 
(Mr) 


d 

(kpc) 


(K) 


iV(CO) 
(cm' 2 ) 


i 

n 


(AU) 


-Rraax 

(AU) 


M(CO) 

(M ) 


xl 


16164-5046 
11097-6102 
18006-2422 
08576-4334 


16164nr3636 
11097nrl218 
18006nr766 
08576nr292 


30 
11 
11 
6 


3.6 
2.8 
1.9 
0.7 


4480 
2710 
1800 
1660 


3.9 x 10 20 
8.9 x 10 20 
1.0 x 10 21 
3.9 x 10 21 


30t^ 

io+i° 

60+ 20 

DU -30 

27l 2 4 


O 1+18.2 

°- 1 -2.9 

6+ 01 

u -°-0.4 
n 9+0-24 
u - z -0.06 

0.2^:1 


+123.5 
°- z -3.0 

q q + 13.5 

°-°-3.1 

1 Q+2-5 
J - ! '-1.6 
o fi +70 


2.0 xlO -9 
2.6 xlO" 7 
6.2 xlO -8 
8.4 xl0~ 7 


0.99 
1.51 
1.72 
2.01 



Table 1. Stellar properties and fit parameters of the first (v — 2 — 0) bandhead. The masses of the objects (column 
3) are based on spectral type estimates, determined from their positi on in the (K, J — K) color-magnitude diagram 
ijBik et al.ll2004|) . The reference of the distances (column 4) are given in lBik et alJ l)2004|) . The errors on the inclination 
reflect the range of inclinations where a good fit was possible. The errors on R m in and R ma x are the corresponding 
values for these inclinations. The errors on T ex and N(CO) are dominated by systematic instead of statistical errors. 



citation temperature T ox , the disk inner and outer radius 
i? m i n and -R max (Table corresponding to the best fit- 
ting model for each object (see the bottom spectra in each 
panel of Figs. ^ and The stellar mass and the inclina- 
tion are degenerate parameters. The parameters are best 
constrained by the first bandhead which is not contam- 
inated by the P-branch of the other bandheads nor by 
Hydrogen Pfund lines. For 11097nrl218, we were able to 
fit the four bandheads (Fig. [21 which shows that T cx is 
lower for higher v (see next section). 

4. Results and Discussion 

We have observed the CO bandheads in a sample of 
15 massive YSO s and five show these lines in emission 
ijBik et alj ]"2004). In this Letter, four of these objects are 
modeled and are well described by optically thin CO gas 
located in a keplerian disk. From the quality of the fits, it 
is clear that additional emission from a disk-wind is not 
required. 



Our results show that the CO gas is located within 
a few AU of the central star (Table . The maximum 
distance where the CO bandheads are emitted depends on 
the inclination, which is not a well constrained parameter, 
but even for a wide range of possible inclinations the hot 
CO molecules are located within a few AU from the star. 

In the absence of extinction by dust grains, the CO 
molecules should be photodissociated by the stellar ul- 
traviolet photons. However, the derived column densities 
(10 20 - 10 21 cm~ 2 ) are well above the required value 
for the CO molecules to self -shield (iV(CO) ~ 10 15 cm- 2 , 
Ivan Dishoeck fc Blacklfl988h . Moreover, gas phase chem- 
ical models show that CO molecules can rapidly form in 
the gas phase to compensate for their destruction. 

The CO column densities in the four objects only dif- 
fer by at most a factor of 10. This reflects the narrow 
parameter range in density and temperature at which the 
CO first-overtone bandheads are emitted. A steep density 
and/or temperature gradient can result in a very narrow 
region in the disk where the CO first-overtone is emit- 
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ted. This is reflected in the value of R m in and R ma x in 
e.g. 16164nr3636, where only a narrow ring in the possibly 
much more extended disk is producing the CO bandheads. 

The derived excitation temperatures range from 1700 
to 4500 K. This temperature may not be equal to the 
gas kinetic temperature for the high-v transitions since 
the critical density is increasing with v (rin > 10 10 cm -3 ). 
In addition, the fit to all four bandheads shows that T cx 
is lower for higher v transitions. Therefore, the popula- 
tion of the high v levels is likely sub-thermal. The exci- 
tation temperature shows a positive trend with increas- 
ing effective temperature of the central star (Table QJ. 
One possibility is that the diluted near-infrared radia- 
tion field pumps the higher u-levels. Likewise, UV fluo- 
rescence will advantageously populate the high w-levels. 
Alternatively, the earlier type stars have higher luminosity 
which results in an increased heating by photo-electrons. 
The later explanation is favored because both radiation 
pumping mechanisms (UV and near-IR) will preferably 
populate the higher rather than the lower i;- levels, which 
is not corroborated by the data. However, only a detailed 
model that includes all the CO levels and possible pump- 
ing mechanisms will give a definitive answer. 

Two scenarios can explain the presence of an inner disk 
at the stage of massive star formation where the circum- 
stellar environment is being dispersed by the stellar wind 
and radiation. 

According to the first scenario, the disk material may 
originate from the rapidly rotating star itself. However, 
only the wind of a B[e] supergiant can produce a disk 
with sufficiently high densities to emit CO bandheads 
l|McGregor et alll988tlkraus fe Lamersl2003)h The stellar 
charac teristics exclude the possibility that they are super- 
giants l|Bik et alJl200^) . 

Alternatively, the inner disk may be the remnant of 
a massive and extended (up to 10 000 AU) accretion 
disk detected in the early phase of star formation (e.g. 
Shepherd et alJ200lUBeltran et alJ2004lGarav fe Lizanol 
19991 IChini et al.ll2004h. I n the model of disk dispersal 
bv lHolleTibach"'et alJ l|l99"lh the photoevaporation of the 
outer part is faster than the dispersion of the inner re- 
gion by stellar winds. This model naturally explains both 
the presence of inner disks and the ab sence of large outer 
disks in most He rbig Be stars (e.g. iNatta et all [2000; 
iFuente et al"ll200dl . 

5. Conclusions & Prospects 

We have modeled the CO first-overtone bandhead emis- 
sions of 4 massive YSOs. The different profiles of the CO 
bandheads are all well explained by a simple keplerian 
disk. Regardless of the gas excitation mechanism, the in- 
formation on the gas kinematics given by the shape of the 
CO bandheads shows that the CO is located close to the 
star (0.1-5 AU). 

The CO overtone emission likely arises in a narrow 
torus where the gas temperature is between 1500 and 
4500 K. This range of temperature reflects the difference in 



effective temperature of the central objects. The decrease 
in excitation temperature for one object from the v = 2-0 
bandhead to the v = 5-3 bandhead indicates that the 
higher u-levels are sub-thermally populated. Excitation 
mechanisms other than LTE, like UV or infrared pump- 
ing and photo-electric heating cannot be excluded by the 
simple model used here. 

The detected disks may be the remnants of massive ac- 
cretion disks. In particular, it is the first evidence of a disk 
around a young 06-08 star, suggesting that these objects 
are recently formed as the disks are not yet dispersed by 
the stellar FUV photons and stellar wind. 

Spatially resolved observation of CO bandhead emis- 
sions with Amber at the VLT interferometer will provide 
further constraints on the size and the geometry of the 
CO emitting region. 
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